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Abstract 

Results are presented for a simulation carried out to test the precision with which a detector design (HERON) based on a 
superfluid helium target material should be able to measure the solar pp and 'Be fluxes. It is found that precisions of ±1.68% 
and ±2.97% for pp and 'Be fluxes, respectively, should be achievable in a 5-year data sample. The physics motivation to aim for 
these precisions is outlined as are the detector design, the methods used in the simulation and sensitivity to solar orbit eccentricity. 



1. Introduction 

This paper presents an analysis of the projected capa- 
bility of a detector design, HERON [ 1 ] , based on a target 
material of superfluid helium to make a precise mea- 
surement of both the pp and 'Be solar neutrino fluxes 
(Opp and <I>Be, resp.) in a single, real-time experiment. 
The detection reaction used would be the elastic scat- 
tering of neutrinos by electrons (ES). In addition to the 
novel use of helium, the detector also includes the novel 
application of a coded aperture [2,3] technique for ac- 
curate measurement of the location and recoil energy of 
each elastic scattering event and to aid in background 
discrimination. According to models of the Sun [4], the 
neutrinos from the pp-I and pp-II branches of the fusion 
chain (known as the pp and 'Be neutrinos, respectively) 
are, when taken together, > 98% of the neutrino flux 
and are associated with the reactions producing a sim- 
ilar fraction of the solar energy. At the present writing 
there have been no real-time experiments to measure 
the flux and spectra of Opp but recently the Borexino 
collaboration has made the first real-time spectral de- 
tection of <I>Be(861 keV) [5]. As we explain in Sec. 2 
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there are several important physics issues which can be 
addressed if a detector can be constructed to measure 
both these fluxes and spectra with sufficient precision. 
For the pp neutrinos (< 420 keV) there does not yet ex- 
ist any detector with demonstrated feasibility to mea- 
sure either their flux or spectra; however, there are a 
number of efforts [1,6] which aim to do so. 

Section 2 of this paper discusses the physics goals 
motivating the HERON detector. Sec. 3 briefly presents 
the requirements for and description of the HERON 
detector design. Sec. 4 and 5 provide the details of the 
HERON capability analysis. In an Appendix we discuss 
an application to measuring the solar orbit eccentricity. 

2. Physics Goals 

A principal goal of HERON would be to make an ac- 
curate measurement of the luminosity of the Sun using 
precise measurements of active neutrino fluxes. An ex- 
periment capable of measuring both <t>pp and <I>Be suf- 
ficiently well for an accurate luminosity measurement 
can also address several other interesting topics. These 
include testing for the relative rates of the ^He(^He, 
2p)^He and ^He^^He, Y)'Be reactions which terminate 
the pp-I and pp-II branches in the Sun and also for test- 
ing the MSW effect (after Mikheyev, Smirnov, Wolfen- 
stein [7]) in the LMA (large mixing angle) solution 
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to the "solar neutrino problem". Additionally, if new 
measurements of Opp and <I>Be are successful at the 
few-percent level then, via new, luminosity-constrained 
global fits to all neutrino data, some modest improve- 
ment can be made in the knowledge of tan^ 9i2, sin^ 613 
and limits on sterile neutrinos. Lastly, since a real-time 
low-energy solar neutrino experiment opens a new win- 
dow in neutrino physics, the possibility of surprises in 
the physics should not be discounted. 

2.1. Solar Luminosity. Why measure the solar 
luminosity by neutrinos? 

The radiant photon energy reaching the Earth from 
the Sun (the irradiance, /g) is believed to result from 
the nuclear fusion reactions of light elements. The en- 
ergy released in each of the chains producing neutri- 
nos is well known from laboratory experiments. Con- 
sequently, if the flux of the associated neutrinos can be 
determined then the photon irradiance and luminosity 
can be inferred from those fluxes. 

This is usually formulated [8] as; 
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= ^(a,<I>,) (1) 



where is the total solar luminosity in photons, 
(A.U.) is the average Earth-Sun distance, Ie is the 
mean irradiance determined by Earth-orbit satellites to 
be 1358.8Wm-2(8.482 X lO^MeVcm-^g-i) with a 
systematic uncertainty of about 0.4% [8,9,10]. The a,s 
are the coefficients giving the energy provided by and 
associated with the i-th neutrino flux <I>, and ^©(v) is 
the total solar luminosity inferred from the neutrino 
fluxes. 

If the Sun operates as we presently think it does 
[ 1 0, 1 1 , 1 2] then the ratio L© (v ) /L© should be unity. Sig- 
nificant departure from that expectation would signal 
the presence of different sources of energy within the 
Sun [13,14,15]. Another important point is that, from 
the reaction positions in the solar interior, the energy 
carried by the photons and by the neutrinos reaches 
the Earth with a huge separation in arrival times. The 
neutrinos arrive directly in 8 minutes while the thermal 
photon energy arrives from the solar plasma after ap- 
proximately 40 000 years [16]. Consequently, finding a 
disagreement between L© and L©(v) would have sig- 
nificant implications for environmental consequences in 
the long term. 

Because the sum of pp (91.5%) and ^Be (7.4%) neu- 
trinos are expected to be associated with > 98% of the 
total flux, it follows that a precision measurement of 



<t>pp, either alone or together with <I>Be, will provide the 
major direct test of L©,(v)/L©. Currently this compari- 
son is only known to about 25% [17,18]. 

There are additional reasons to make a more precise 
determination of L©(v)/L©: the fact that the average 
photon irradiance is very well measured [9,19,20] has 
previously led to its use as a constraint in global anal- 
yses of solar neutrino experimental data. Used first 
as a demonstration of possible flavor oscillations of 
neutrinos [21] and more recently, as additional and 
more precise solar and reactor neutrino experimental 
results have become available, as a powerful constraint 
to aid in establishing best present knowledge of solar 
neutrino mass-mixing parameters and individual fluxes 
[17,22,23,24,25,26]. 

At the same time, the Standard Solar Models (SSM) 
[10, 11, 12] have continued to make significant improve- 
ments so that quite precise predictions for the fluxes 
have been made. For example, Opp is predicted to ±1% 
and <t>Be to ±9.5% [llj3- When these predictions are 
compared against the fluxes found from global fits to 
all of the existing solar and reactor data [28] the levels 
of agreement differ significantly depending on whether 
the photon luminosity is used as a constraint or not. 
For example, the ratios of global-fit fluxes to SSM 
predictions are (at la): for pp, 1.01 ±0.02 and for 
^Be 1.03^j Q3 with the luminosity constraint but are 
1.38^Q25 0.13^Q^3, respectively, without the con- 
straint and leads to the poor knowledge of L©(v)/L© 
noted above [17,24,25,26,29]. 

The question of how precisely direct measurements, 
of either Opp or <I>Be must be made in future experiments 
has been cogently addressed in an important paper by 
Bahcall and Pena-Garay [17]. Related and more recent 
considerations also have been made by others [11,18]. 
The level of precision required depends strongly upon 
the specific physics questions to be addressed. In all 
cases, the demands on experimental techniques are se- 
vere. For example, the authors of Ref. [17] carried out 
simulations of global analyses utilizing all present data 
plus inclusion of potential future pp and ^Be experi- 
ments with assumed capability of precisions ranging 
from 1% to 30% (at la). They find that a ''Be result 
of ±5% could improve knowledge of L©(v)/L© from 
1.4^Q3 to 1.07±0.13. Increased precision on ^Be alone 
would not yield further improvement; while a ± 1 % on 
pp could achieve a remarkable 0.99 ±0.02 on the lumi- 



' Note 1: 8.5% of the 9.5% is contributed from the expeiimental 
uncertainty in 534 nuclear cross-section factor; however, new data 
from the LUNA collaboration [27] suggests this contribution may 
be reduced to 2.5%. Private communication C. Pefia-Garay. 
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nosity comparison by neutrinos. The authors note that 
a resuh of this accuracy would be "a truly fundamental 
contribution to our knowledge of stellar energy gener- 
ation and place a ±2% bound on all sources of energy 
other than low energy fusion of light elements (i.e., pp 
and CNO chains)". The HERON detector, as shown in 
the present paper, is intended to be capable of reach- 
ing precisions on both fluxes commensurate with these 
goals. 

2.2. The pp-I vs pp-II and LMA-MSW 

The relative magnitude of <t>pp versus <t>Be is a par- 
ticularly relevant parameter, on the one hand, bearing 
on the accuracy of the SSM and, on the other, as evi- 
dence for the MSW effect in the LMA. It is valuable 
to have an experiment which measures both Opp and 
4>Be since several systematic errors tend to cancel in the 
ratio <t'Be/*I'pp- In the SSM there is a very strong anti- 
correlation between the two fluxes with a coefficient of 
-0.79 to -0.81 [11]. If the 3HeC^He,Y)^Be reaction 
of the pp-II chain were the only terminating branch, 
only one pp and one ^Be neutrino would be produced 
in each cycle. Otherwise there would be two pp and no 
^Be neutrinos if ^He(^He,2p)^He of the pp-I branch 
were the terminating reaction of the full fusion cycle. 
What the actual relative reaction rates are depends on 
several not yet accurately known details within the Sun 
such as elemental abundances, temperatures and den- 
sity. (Present versions of the SSM predict a ratio of the 
two reaction rates as 0.174 [24] which implies, prior 
to oscillation, a value of 0.080 for "t>Be/"t'pp ) An inde- 
pendent and precise measurement of these relative rates 
would be an important contribution to the understand- 
ing of stellar processes and would permit a refinement 
of the use of the SSM in global analyses of neutrino 
data. The physics of the MSW effect [7] is embodied in 
the flavor-dependent interaction differences for neutri- 
nos propagating in matter as opposed to vacuum. Due 
to the differences in neutrino energies and solar density 
at their production points the pp, ^Be and ^B neutrinos 
are expected to have quite different survival probabil- 
ities. The LMA-MSW solution specifies what this en- 
ergy dependence must be. 

The oscillations of the much higher energy ^B neu- 
trino should be strongly suppressed by matter domi- 
nance and the pp neutrinos much less since they should 
be vacuum-dominated. The ^Be and pep neutrinos, hav- 
ing energy intermediate to pp and ^B, are in the cru- 
cial energy region where the transition between matter 
dominance and vacuum oscillations is to be expected. 



The B flux is now very well measured by the Super 
Kamiokande (SK) and Sudbury (SNO) experiments [in 
Ref. [28], see e.g. Fukuda et al. and Ahmed et al.]; how- 
ever, direct experimental evidence for this MSW tran- 
sition is still lacking and could be provided by an ex- 
periment such as HERON. 

2.3. Other areas of interest related to neutrino 
properties 

The Opp and <I>Be measured in ES are, by necessity, 
the fluxes of active neutrinos. Consequently, in the mea- 
surement of Lq; (v)/L0 an alternative interpretation of a 
result consistent, within errors, to unity can be taken as 
the establishment of a limit on the presence of sterile 
neutrinos. 

Due largely to the lack of precision experiments 
on these two major low-energy fluxes, there has been 
some leeway in the recent analyses of the solar and 
reactor data which allows for consideration of sev- 
eral well-motivated proposals for "new physics" in the 
neutrino sector. Among these are possibilities for non- 
standard neutrino interactions with their environments 
(NSI) [30]. Within this class of models the additional 
effects to be expected are strongly constrained by exist- 
ing data; however, in some cases they should be most 
pronounced in the energy dependence of the fluxes of 
solar neutrinos < IMeV. In these cases the effect is 
qualitatively similar, but differs quantitatively, from 
that to be expected from the LMA-MSW transition 
in the matter- to vacuum-dominated neutrino energy 
regions. Two examples of this class, are models with 
flavor-non-conserving neutral currents [31] or with 
mass-varying-neutrinos ("MaVaNs") [32]; the latter 
inspired by possible insights into understanding "dark 
energy". Confirmed evidence found for NSI would 
place our knowledge for the mass-mixing parameters 
in doubt; alternatively, such precision measurements of 
the matter-vacuum transition region would also serve 
to establish new limits on the existence of NSI. 

New measurements of low energy solar fluxes can 
play only a rather limited role in improving on present 
knowledge of tan^0i2 and sin^9i3 limits. To be useful, 
the new data would need to be folded into a compre- 
hensive analysis with all the other data (solar, atmo- 
spheric and reactor) which from present data have al- 
ready established impressive errors on the parameters 
[18,23,26,33]. The potential for improvement in these 
parameters by low energy ES or CC experiments has 
been subjected to a detailed study by the authors of 
Ref. [17] who conclude that without new physics even 
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± 1 % on <I>Be would make a negligible improvement on 
tan^ 9i2 and a <J>pp result < ±3% is required to improve 
the error by more than 15%. The authors' conclusions 
on any improvements to be expected on sin^0i3 are 
similar. 

3. Detector 

3.1. Requirements 

There are stringent requirements placed on any de- 
tector designed to achieve these goals. In order to be 
sensitive to all active neutrino flavors, the detection re- 
action is that of the elastic scattering from atomic elec- 
trons (ES) in the target: Ve,;/,x + e^ ^ Ve./^,.x + e^ . The 
ES event signature is the occurrence of only a single, 
low energy recoiling electron in the detector medium. 
The recoil spectra are continuous from zero with the 
pp decreasing monotonically to a maximum energy of 
26 Ike V while the ^Be spectrum is nearly flat up to 
664 keV maximum (see points labeled input in Fig. 4). 
These recoil spectra place a premium on achieving as 
low an energy threshold as possible. The most danger- 
ous backgrounds are those which can be created by 
the appearance of electrons from Compton recoils of 
gamma rays or radioactive decays within the medium 
of the target. These backgrounds need to be mitigated 
by a combination of event signatures, use of low activ- 
ity materials, depth and shielding. In order to achieve 
the desired few-percent precision, very large statistics 
event samples are required and systematic errors arising 
from effects due to analysis cuts (e.g., fiducial volume, 
thresholds, pp and ^Be event separation, calibrations) 
must be minimized. Although these details emphasize 
the challenges to be faced in constructing an operat- 
ing detector, there are two mitigating factors: the ES 
cross-section is very precisely known due to experiment 
and electro-weak theory and the expected fluxes are 
high (4>pp = 5.99 X lO"* and <I>Be = 4.84 x 10''cm-2s-i 
[11]). Typically, this will result in an ES event rate of ~ 
2 events • tonne-day^' which implies that only a modest- 
size fiducial volume (e.g., ~ lOtonne) is needed for a 
high statistics experiment. 

3.2. HERON configuration and detection processes. 

The detector design is shown in (Fig. 1) and its de- 
sign is discussed in more detail in Ref. [1,34,35] . The 
general approach of HERON to these issues is as fol- 
lows. The target material chosen is "^He in the superfluid 
state (density 0.145 g/cc) which has several beneficial 




Fig. 1. The geometry of the HERON detector design, all dimensions 
ai'e in centimeters. 

properties. Energy deposited in the helium by recoiling 
particles can be detected by one or all of three pro- 
cesses: scintillation, phonons/rotons or collecting the 
recoil electron trapped in a bubble. Helium has no long- 
lived isotopes but more importantly it can be made ab- 
solutely free of all other atomic species. At superfluid 
temperatures it is self-cleaning of impurities due to their 
high mobility and favorable energy minimum at the con- 
tainer walls. Even particulate matter quickly attaches to 
the container walls at our operating temperature of ^ 
30-50 mK. Since the bulk helium volume will be free 
of background sources, the concern is to counter ra- 
diation entering from the cryostat (27.3 tonne copper) 
and its environment. In a separate study environmen- 
tal sources were modeled with the detector cryostat ac- 
tively shielded externally by 3.5m of water at a rock 
overburden of 4500 meter-water-equiv. (m.w.e.). He- 
lium is virtually immune to creation of long-lived cos- 
mogenic muon activity, capture or decay in it; muons 
(4m-2day"' at 4500 m.w.e.) are vetoed externally and 
internally in any case and greater depth is possible. 
This was sufficient shielding against environmental neu- 
trons, cosmic muons and gammas that, for purposes of 
the analysis simulation under consideration here, they 
would be negligible relative to sources from the cryostat 
and other detector parts. As a consequence the back- 
ground issue reduces to controlling conversions of gam- 
mas entering the helium volume of 2 1 .6 and 8 tonnes to- 
tal and fiducial, respectively. However, helium does not 
have good self-shielding properties and this is partially 
compensated for by lining the cryostat with a moderator 
of solid nitrogen enclosed in acrylic cells (5.6T acryUc 
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and 104T solid nitrogen, density 1.03g/cc). The func- 
tion of the moderator is to absorb or degrade in energy 
by Compton scattering the entering gamma rays which 
originate 97% and 3% from the cryostat and moderator, 
respectively. The flux of gammas entering the helium is 
dominated by low energy (dominantly < 1 MeV) cos- 
mogenic activity in the copper with the remainder from 
U, Th and other activity in the copper, nitrogen, acrylic 
and other parts. The background gamma conversions in 
the helium cannot be fully eliminated but are amenable 
to the development of a distinguishing signature which 
aids in their separation from signal; the nature of this 
separation is detailed in Sec. 4 and 5. 

We have chosen to use as the basic elements of de- 
tection for both signal and background the collection 
of scintillation light and also the recoil electron trapped 
in its bubble. We have carried out studies of these pro- 
cesses using prototype calorimeter sensors/detectors de- 
veloped for this project and suitable for both signal 
types; the processes and results are discussed briefly 
below and in more detail in the references cited. 

Excited or ionized He atoms along the electron path 
quickly form dimers in the liquid. The radiative decays 
to the ground states of these singlet and triplet dimers 
emit photons in the ultraviolet. The scintillation light is 
in a narrow band centered at 16 eV and results from 
the decay of the singlet dimer He2(A'E+) [36]; since 
this energy is lower than the first excited state of He at 
20.6 eV, the liquid is self-transparent. 35% of a recoil 
electron's energy is released 27000 photons/Me V) in 
this singlet dimer mechanism [34]; the energy from the 
long-lived (~ 15 sec) triplet dimer escapes or is colli- 
sion quenched. (An additional 43% of the recoil energy 
is radiated in phonons and rotons which by quantum 
evaporation [37] could in principle also be utilized with 
the calorimetric sensors for discrimination [38]; how- 
ever, we find incorporating into event signatures the de- 
tection of the recoil electron from the drifted bubble a 
much stronger discriminant.) 

When a recoil electron in He has lost most of its en- 
ergy it forms an electron bubble. The electron experi- 
ences a strong, short-range repulsive potential from the 
bound electrons on surrounding atoms due to Pauli ex- 
clusion. This repulsion forms a vacant volume, or bub- 
ble, of 19 A radius in which the electron is confined. 
The bubble forms in about 10 pico-sec with an effec- 
tive displacement mass of ^ 500 He-masses and has a 
hydrodynamic mass of half that [39]; consequently due 
to this difference in masses, under gravity and with no 
electric field, a bubble would experience a buoyant ac- 
celeration of ~ 2g. A uniform drift velocity of the bub- 
ble can be provided and controlled by a combination of 



applied electric field and a very low concentration of 
^He; for example, at 40 mK, 30ppm ^He and a field of 
300 V/m provides a drift velocity of 17.5 m/s. In a worst 
case example of 5 m (maximum depth), the 3x10^ col- 
lisions induced due to the large cross-section of ^He for 
scattering a bubble leads to an uncertainty in the transit 
time of 16^s and hence to a depth error of < 1 mm. In 
addition to providing a "drag" force the ^He also aids 
in extracting the electron efficiently through the free 
surface of the liquid by vortex attachment (Surko and 
Reif [40]). Two grids on either side of the liquid sur- 
face (not shown in Fig. 1) provide the drift and extrac- 
tion fields [43]. (At the normal operating temperature of 
30-50 mK the vapor pressure is sufficiently low that the 
space above the liquid is effectively a vacuum.) The fi- 
nal grid accelerates the electron to ^ 1 keV thus provid- 
ing a large and distinguishable pulse in the calorimeter. 

For several reasons, photo-multiplier tubes (PMT) are 
not suitable event detectors for HERON, among them: 
the high radioactivity of PMT's, poor He self-shielding, 
lack of transparency of moderator and the desire to de- 
tect the drifted electrons. As mentioned, both scintil- 
lation and drifted electrons are detected on the same 
calorimetric devices. Each device constitutes a pixel in a 
geometric array (a coded aperture) and consists of a thin 
wafer of silicon or sapphire to which is attached a high 
sensitivity metallic magnetic calorimeter (MMC) read 
out with a SQUID sensor [41]. For astrophysical x-ray 
application, versions of MMC have been constructed 
with ~ 3eV resolution [42]. Projecting from measure- 
ments on wafer prototypes of smaller heat capacity to 
ones of the HERON size, ~ lOeV resolution is to be 
expected. In a full simulation of the response of this 
large wafer (100 x 100 x 0.4 mm, in this example) with 
16 eV photons it is found that single photons should be 
detectable at wafer temperature of 40 mK producing a 
pulse of 5 ms rise- and 100 ms fall-time [44]. This per- 
formance capability is assumed in the context of the 
analysis of Sec. 4 and 5. 

For each neutrino event we must reconstruct its po- 
sition within the He and also its recoil electron energy. 
In addition, we must develop event signatures which 
aid in separating signal from background events. The 
maximum track length expected for a neutrino event is 
^ 2 cm so that on the scale of the total helium volume 
(149m^; 21. 6T) neutrino signal events are effectively 
point sources of scintillation light. At the dominantly 
low energies of the gamma-ray background events the 
conversions in He are overwhelmingly (95%) Comp- 
ton scatters. 90% of these conversions are multiple de- 
positions distributed over an average length of more 
than 50 cm in the He. Consequently, the scintillation 
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from background arises from a distributed, rather than 
a point, source; additionally the event most often con- 
tains multiple, un-recombined electron recoil bubbles. 
These latter two features constitute the primary back- 
ground signature. Subsequently, differences among the 
spectral and spatial distributions of the signal and resid- 
ual background events facilitate their final separation. In 
order to create these signatures and to enable the neces- 
sary cuts on data samples the 2400 wafer calorimeters 
are arranged into two planes in the vacuum space above 
the liquid; the resulting array constitutes a coded aper- 
ture and provides the ability for both spatial and energy 
reconstruction. 



4. Nature and use of HERON's coded aperture. 




4.1. General principles. 



Fig. 2. URA pattern, 40 x 40 in dimensions and 19 x 17 in periodicity. 
Dark and light areas represent opaque and open regions, resp. In 
the HERON mask, the smallest squares represents single wafers. 



The concept of coded aperture arrays [2,3] has been 
a well established one with arrays being widely used in 
x-ray astronomy . The role of the array is to accurately 
determine the direction of incoming photons from a re- 
mote x-ray point source. A coded aperture array con- 
sists of two parts, an imaging plane and, separated by 
a fixed distance, a mask plane. In the x-ray application 
the image plane consists of a set of active sensors (or 
pixels) while the mask plane is opaque with a pattern 
of cut-out apertures. For far-distant sources, nearly par- 
allel rays enter the array and some of them are blocked 
by the opaque portions of the mask; thus the image im- 
posed on the sensor plane will resemble the pattern of 
the apertures in the mask effectively leaving a "shadow" 
(See Fig. 2). In principle, to determine the direction of 
the light source it is simply a matter of comparing the 
shadow pattern to that of the mask itself. In practice, 
elegant techniques have been developed for design of 
mask patterns and image deconvolution taking into ac- 
count side-lobe effects as well as intrinsic and statistical 
noise. Various classes of mask patterns have been em- 
ployed in the x-ray field ranging from random apertures 
to strictly repeating and regular patterns. The choice of 
one over another depends upon experimental considera- 
tions such as strength of photon flux, resolution needed 
and sidelobe tolerance. 

In the HERON application there are important dif- 
ferences: a) a non-distant point source location is to be 
determined within a limited volume in 3-D, b) the event 
energy must be measured, c) for reasons (a) and (b) both 
planes will consist of active wafer pixels and d) some 
discrimination is needed between point and non-point 
sources. The HERON coded aperture array is arranged 



with the mask plane 1 cm above the liquid surface and 
the image plane Im above the mask. The pixels are 
arranged in 40 x 40 square arrays (460 x 460 cm each) 
with 1600 wafer calorimeters in the image plane and 
800 in the mask. The pattern of apertures in the mask 
is that of a uniformly redundant array (URA) [3]. This 
mask pattern was chosen because of its low intrinsic 
noise and high transparency (50%). (In the notation of 
the URA it has a (17,19) grid spacing [3] as shown in 
Fig. 2.) The nature of the URA and HERON physical 
properties and goals constrain the choice of wafer pixel 
size. The He fiducial volume can be chosen and varied 
during physics analysis but the total mass of He is con- 
tained within a cylindrical volume of R = 330cm and 
435 cm height. In our application, the transverse dimen- 
sions of the array should be at least commensurate with 
those of the fiducial volume containing the sources and, 
although smaller pixel sizes imply finer spatial resolu- 
tion, ultimately photon statistics dominate (typically a 
few hundred photons for many neutrino events, due to 
energy and solid angle effects). For the HERON geom- 
etry a choice of 11.5 x 11.5 cm pixels in a URA gives 
resolution adequate to the physics goals without unnec- 
essarily increasing the complexity or noise and consis- 
tent with the desired single photon performance for a 
wafer of this size. 



4.2. Loglikelihood method for event reconstruction. 

The deconvolution techniques used for typical x-ray 
applications are not applicable for our 3-D application; 
additionally, they are not easily amenable to developing 
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a background (distributed source) signature. Instead we 
have adopted a likelihood approach along with a search 
algorithm for the most probable position in 3-D space. 

This approach treats each event in a sample contain- 
ing both signal and background as if it were a single 
point source. With that assumption, it finds from the 
observed photon hit pattern the most probable values 
of its spatial location and total energy. Although no at- 
tempt is made on an event-by-event basis to distinguish 
signal from background, distributions of the likelihood 
parameter's logarithm can be useful in separating sig- 
nal and background as we show in Sec. 5. Similarly, 
the effective point-like positions and effective energy 
distributions of the background events are used. 

Operationally, the algorithm initiates with a test-point 
location in the volume and the probability of this test 
electron to produce the observed photon hit pattern 
is calculated. The test photon distribution is taken as 
isotropic with straight-line propagation; the probability 
of hitting the i-th wafer is then proportional to the solid 
angle (co,) subtended by the wafer from the photon cur- 
rent starting point. After a systematic search of points 
throughout the available volume, the location in space 
found to have the highest probability is taken as the fi- 
nal position. 

If Q. is the solid angle subtended by all m wafers 
in both planes and is the number of photon hits in 
the pattern then we can define a quantity A'tot = 4-nN/Q. 
which is the total number inferred for the test point. 
Then the probability of an event located at x producing 
the recorded photon pattern is evaluated as: 



^ 1 \ jVtot m 

"(«)=! 5;) n 



CO;(x)"' 



(2) 



i=0 



where «, is the number hitting i-th wafer and «o = 
Mot — N and for computational convenience we use the 
logarithm (loglikelihood); 



L (x) = InP(x) = -Motln47i: + £ In 



co,-(x)"' 



(3) 



and select as the final position the one with the largest 
(least negative) L . The final energy estimate is scaled 
from the solid angle subtended from the test point. The 
process converges rapidly guided by a set of empirically 
established criteria for avoiding subsidiary maxima and 
reaching a stable solutionis 



^ One might expect tlie determination of tlie effective solid angle 
and efficiency to depend on the optical properties of the UV radiation 
for He and the wafers. Effects resulting from: the width of the 
singlet dimer spectrum, internal reflection and refraction at the He 
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5x lO-'^g/g 


4x lO-'g/g 





Table 1 

Assumed levels of residual activity in major detector components. 
4.3. Event simulation and reconstruction 

A test of the reconstruction ability of the coded aper- 
ture approach has been done in a way which examines 
a full range of variable correlations. We have gener- 
ated samples of pp, ^Be (5 x 10^, each) and background 
events (2.5 x 10^) as they would appear in the configu- 
ration of HERON described. For the neutrinos, the in- 
put recoil electron energy spectra are as shown in Fig. 4 
and the events are distributed uniformly throughout the 
full He volume. The input background sample in the He 
is generated by propagating gamma rays initiating from 
sources within the detector's principal components us- 
ing GEANT3 [47]. The source activities and concentra- 
tions are listed in Table 1. Within the He account was 
taken for bremstrahlung, very low energy Compton re- 
coils and delta-rays. For all input samples, the original 
position and deposited energy for every recoil electron 
was retained for use in comparing to reconstructed val- 
ues. The events were then reconstructed as described in 
Sec. 4.2. 

Signal reconstruction. The principal interest in this 
section concerns the reconstruction results for the true 
point sources — the pp and ^Be events. 

Some examples of spatial and energy resolutions (1- 
o) for the reconstructed pp and ^Be samples are shown 
in Fig. 4, 5, and 6. Rxy- (calculated from the fitted x and 
y coordinates) is the horizontal position in the He cylin- 
drical volume and E is the total deposited energy. The 
7-coordinate (depth) is a factor ~ 2 x poorer when re- 
constructed by coded aperture; however, since the bub- 
ble drift time gives better than 1 mm resolution the depth 



surface, Rayleigh scattering, multiple reflections between He and 
wafers or among wafers themselves can affect the value of the 
effective soUd angle. However, these effects are small due to the 
low refractive index of He (1.045), the self-transparency and very 
long Rayleigh mean-free-path (200 m [45]), the narrowness of the 
UV spectrum and the low reflectivity of the wafers at the incident 
angles involved [46]. These effects can be taken into account by 
including ray-tracing techniques and we have tested them in our 
analysis [35]; however, it is more computationally demanding and 
results in no essential differences on flux as we demonstrate in our 
final error table (Table 5). Consequently, for simplicity of analysis 
and discussion we ignore these effects in this paper 
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position by timing is used instead. While the resuhs of 
the reconstruction are quite satisfactory, several points 
should be noted about the figures and the nature of the 
results. The resolution dependence on photon statistics 
is common to both energy and position and is dominant 
at the lowest energies, but it can also arise at higher en- 
ergies as a contribution due to diminished photon solid 
angle from the extremities of the detector and mask ge- 
ometry. As the two inserts for Fig. 5 and 6 show, the 
maximum value of the distributions indicate excellent 
agreement between input and output values, the shapes 
near the maxima are symmetric and the FWHM look 
reasonable. However, such figures are too Umited to ad- 
equately illustrate the effects of correlations nor of any 
asymmetries in the tails. The correlations between spa- 
tial and energy variables are of major concern as they 
enter into flux accuracy when making the high level 
analysis cuts for fiducial volume and threshold. For our 
main purpose, the precise neutrino flux determination, 
it is most direct and useful to examine comparisons of 
the complete set of distributions to be used during the 
flux separation of the three channels of events. This is 
the case because if our model of the detector is correct 
the full generation and reconstruction of these samples 
naturally includes all of the effects of the correlations 
and tails and permits estimation of their contributing 
errors. Sec. 5 examines the question of cuts to reduce 
background and to make flux separations. 

Background reconstruction. Any event trigger is ac- 
cepted with a scintillation pulse above hardware thresh- 
old. An important point in the subsequent treatment 
of the events is dependent upon the detection of the 
electron(s). The difference in electron bubble topology 
(multiplicity) between background and signal is an ad- 
ditional key factor in reducing the raw background rate 
for analysis. 

The nature and magnitude of the sources used in our 
simulation were listed in Table 1 . Any event giving one 
or more gamma conversions anywhere in the entire he- 
lium volume is recorded and reconstructed by the coded 
aperture method of Sec. 4. The simulation result im- 
plies a total daily rate in the full He volume of nearly 
4x10^ events; a factor ~ 70 greater than signal. Several 
properties of the background depositions are useful to 
establish prior to making high level cuts on the samples 
for flux separation. These properties (multiplicity, event 
loglikelihood, energy and spatial distributions) were de- 
termined from the sample of 2.5 x 10* simulated events. 

Due to the large size of the He volume and the low 
energies of the gammas entering through the moderator, 
the detector is essentially hermetic and in most cases 
fully contains the event. Because of the He density and 




E (keV) 



Fig. 3. Energy spectra of simulated background events. "Input" 
contains all events before reconstruction; "Fit" are the reconstructed 
energy of these events; "1+2" are the events with no more than 2 
drifted electrons; "1+2 & cuts" are the "1+2" with "std" fiducial 
and 45keV energy threshold cuts. 

gamma mean free path a typical event, prior to cuts, 
consists of widely distributed Compton recoils with a 
mean multiplicity of 12.9 electrons and a mean spa- 
tial separation of 57.6cm. The uncut energy deposition 
spectrum both before and after reconstruction is shown 
in Fig. 3. It should be emphasized that it is not important 
in our method to extract the true energy of each back- 
ground event, which in any case cannot be done with 
a coded aperture for events of these topologies since 
there is no true point-origin of scintillation. In addition 
to providing a non-pointlike i -value, an effective en- 
ergy is obtained and this accounts for the smoothing 
seen in the reconstruction of the raw input spectrum for 
gamma conversions in Fig. 3. 

In order to use the number of detected drifted elec- 
trons effectively there must be a high efficiency for col- 
lecting each one. Full 100% efficiency is not possible 
due to ion re-combination. As noted, for drift control 
purposes we envision a uniform field of modest strength 
(300V/m). With this field and in our geometry very high 
efficiency (~ 85%) can be expected for any electron 
with >5keV [43]. 

As a first step in reducing the background sample for 
further analysis we impose an event selection criteria of 
multiplicity < 3 detected electrons; each with an energy 
> 5keV. This reduces the background by a factor 0.38 
without affecting the neutrino signal (possible brem- 
strahlungs from neutrino events remain included by this 
cut). See Fig. 3 for the resulting effect of this multi- 
plicity cut on the spectrum. Since the goal is to extract 
pp and Be fluxes separately, we choose not to make an 
event-by-event division of neutrino events from back- 
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Fig. 4. Comparison between input neutrino spectra and their recon- 
structed counterparts. Each spectrum contains ~ 5 x lO' simulated 
events and is normalized. 
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Fig. 5. Resolution of reconstructed energy for neutrino events with 
respect to the energy of input events (asterisks: pp; diamonds: Be). 
The inset shows the distribution of reconstructed energy of events 
between 100 and 120keV. 

ground with the HkeHhood ratio. Instead we defer to 
making the simultaneous separation of all three chan- 
nels by the method described in Sec. 5. Additional sam- 
ples of 8 X 10^ background events with detected multi- 
plicity < 3 are also generated. These, together with the 
neutrino samples, are then binned to form PDF's and 
used for further study of high level cuts. 

5. Flux separation 

Fundamental to the method of flux separation used 
here is the construction of accurate probability distribu- 
tion functions (PDF). They present the expected appear- 
ance for four major reconstructed variables {R^, z, en- 
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Fig. 6. Resolution of reconstructed radial position for the neutrino 
events depending on their energies (asterisks: pp; diamonds: Be) , 
The inset figure shows the distribution of reconstructed radii of events 
at different horizontal distances from the center of the detector. 

ergy and L (x)) for events from the two neutrino signal 
channels and the backgrounds. The PDF's are used in 
an extended-maximum-likelihood method [48] applied 
to a simulated experimental random sample including 
events from all three channels. The PDF variables used 
are the radial and vertical positions of the events, their 
energies and the L of Eqn. 3. Our principal goal here 
is to examine the type of systematic errors introduced 
by this approach. 

5.1. Fitting fluxes using likelihood method 

After we have observed and reconstructed neutrino 
and background events with sufficient statistics, we may 
separate the three fluxes (pp, ^Be and background) using 
the extended likelihood method. Given certain binning, 
the extended (log)likelihood function is defined for each 
individual variable as 



3 3 



(4) 



a=l 



bin i 



a=l 



where a denotes the flux type, pai the probability of flux 
a populating the bin ;, «, the total number of events in 
bin /, and A^a the number of events being of type a. A^oi 
are the parameters that are to be adjusted to maximize 
L<^, with the constraint that their sum equals to the 
observed total number of events. 

The derivation and structure of error for each fitted 
flux A^a by this method is discussed in ([35] Appendix 
A). The eiTor matrix V , defined as Vap = (eaSp), in this 
context evaluates to 
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-N + P" 



where elements of matrices N and P are 
Nap = §apA^a and P„p = 

bin i 



(5) 



(6) 



0.25 



where S^p is the Kronecker delta, with indices a and p 
referring to flux types. 

Note that the elements in matrix P may have very 
large values if some bins have very low counts (i.e. 
«, being small), which will lead to reduced numerical 
stability of P^' and consequently of V. To avoid this, we 
chose to use a variable-width binning where bin edges 
are chosen so that each bin contains roughly the same 
number of events. 

Matrix P describes how much one PDF resembles an- 
other. If for two fluxes a and p, pai and pp, are strongly 
correlated over index /, P will have large non-diagonal 
elements, resulting in large valued elements in P^^ and 
V. In other words, if two PDFs are similar in shape, the 
expected errors of their fitted fluxes will be high. This 
property dictates the variables we can use for flux sepa- 
ration: reconstructed R and z distributions are basically 
identical for both types of neutrinos, thus they are not 
suitable for separating ^Be from pp; while energy and 
loglikelihood of reconstructed events provide satisfac- 
tory PDFs for this purpose. 

5.2. The PDF distributions. 

In generating the PDF's, the input pp and 'Be neu- 
trino oscillated fluxes to be expected at Earth are based 
on the Bahcall-Pinsonneault solar fluxes [11,12], the 
mass and mixing parameters from the global fit of Fogli 
et al [23] and assumed two-flavor mixing. 5 x 10^ re- 
coil electron events were generated for each of pp and 
'Be and each contained the appropriate mix of active 
neutrino flavors (Ve, V^x). (The resulting daily event rate 
above 45 keV is 1.7 pp and 0.6 'Be per tonne.) The 
points labeled "input" in Fig. 4 show the spectra in 
HERON to be reconstructed. For background 2 x 10^ 
gamma rays due to the sources listed in Table 1 were 
generated and propagated via GEANT through the de- 
tector. Of these, a sub-sample of 8 x 10^ with less than 
three conversion electrons in the He each with ener- 
gies greater than 5keV were used in creating the back- 
ground PDFs. This sub-sample choice resulted from the 
detailed study of the properties of a full sample of con- 
versions as outlined in Section 4.3. Shown in Figures 
7-10 are samples of the four PDF's prior to any fiducial 
or threshold cuts. 
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Fig. 7. Normalized PDF for reconstructed energy spectra of simulated 
events before high level cuts. These include 8 x lO' background 
events, and > 5 x 10^ events in each neutrino flux. 
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Fig. 8. Normalized PDF distribution of fitted loglikelihood of sim- 
ulated events before high level cuts. 



0.06 

0.05 

S 0.04 
-o 

^ 0.03 

a 

o 

Z 0.02 
0.01 





10 



12 



R (m-) 



Fig. 9. Normalized PDF distribution of horizontal radii fi = \/ + y- 
for the simulated events before high level cuts. 
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Fig. 10. Normalized PDF distribution of depth of simulated events 
before high level cuts. See Fig. 1 for tlie coordinate system. 

5.3. High-level cuts. 

Even after reducing background rate by rejecting 
events with three or more recoil electrons, the back- 
ground conversion rate in the full helium volume would 
still be a factor ^ 25 higher than that of the neutrino 
recoils in the full helium volume, this factor can be 
significantly improved by instituting fiducial volume 
and energy threshold cuts. 

The following considerations were made for choosing 
to impose these high level cuts on fiducial volume and 
threshold with a minimum of bias. Depending on the 
fraction of detector surface area covered by the coded 
aperture, the resolution worsens somewhat at the ex- 
treme edges. In practice the position and energy reso- 
lution will have to be carefully calibrated by temporary 
insertion of sources at measured locations; however, 
cutting in the simulation avoids possible artificially in- 
troduced features in the PDF's due to inaccuracies in 
the reconstruction model. Additionally, it is convenient 
to work with a better signal to noise ratio provided that 
the distinguishing characteristics among the PDF's are 
not compromising the sensitivity. A check on this can 
be seen in Table 2 and 3 where the results are presented 
for analyses with three quite different fiducial volumes 
and two differing thresholds. The contribution to flux 
error is not significantly changed amongst them. Our 
goal has been to seek an optimal set of cutting parame- 
ters that yield the lowest relative systematic uncertainty 
in the fitted neutrino fluxes. 

The various components of errors due to the flux sep- 
aration process and the high-level cuts are discussed in 
the next section. For this purpose we have considered 
a sample size equivalent to 5 years of running in the 
HERON design. This would result, using the "standard" 



fiducial volume (see below), in 2.7 x 10^ pp, 9.3 x 10^ 
Be and 7.5 x lO"* background events. However, in ex- 
amining cut effects we studied six combinations of en- 
ergy threshold cuts and fiducial cuts using the new 
PDF's appropriate to the cut combinations involved. 
The energy thresholds tested are at 30 and 45keV, and 
the fiducial cuts paired with them are referred to as 
"small" (- 177.5 cm < z < 77.5 cm, R„ < 230cm), "std" 
(-187.5cm < z < 87.5cm, Rxy < 250cm), and "big" 
(-208.5cm < z < 108.5cm, Rxy < 270cm), whose co- 
ordinate system is shown in Fig. 1. Interestingly, even 
though the S/N ratios can be very different depending 
on the cuts, the overall errors of fitted neutrino fluxes 
are much less sensitive. 



5.4. Errors in flux separation. 

Besides the statistical error on the total sample there 
are systematic errors in the flux separation process: er- 
ror due to volume and threshold cuts as a result of 
spatial and energy resolution of the reconstruction, er- 
ror from intrinsic uncertainty of the likelihood fitting 
method (arising as discussed in Sec. 5.1) and error due 
to uncertainty of absolute energy scale. We treat the lat- 
ter error separately from the others. 

Error due to volume and threshold cuts stems from 
the inaccuracy of spatial reconstruction of events. Given 
a set of events, one can expect some Na of them to 
be mistakenly accepted, while some other A^^ of them 
to be mistakenly rejected. A'^, and A^,. can be modeled by 
comparing reconstructed events with their correspond- 
ing input values from our simulated event data, and the 
difference A'^, — A^, can be considered as a form of sys- 
tematic bias of fitted flux, which could in principle be 
corrected for but the correction would be very small and 
of the same order as its error Both Na and A^^ are sub- 
ject to statistical fluctuations and such a correction will 
have an uncertainty o = ^/Na +Nr. This error depends 
on both spatial and energy cuts, and its values for some 
of their combinations are listed in Table. 2. 

Errors of separated fluxes inherent to the likelihood 
method, as discussed in Sect. 5.1, are also calculated 
based on simulated neutrino and gamma event data. 
Since these errors depend on the shapes of the PDFs, 
which in turn depend on fiducial and energy cuts, they 
are listed in Table 3 in relation to the cuts also. The 
distributions of fitted energy are used to provide PDFs 
in the table. 

It is worth noting that we are using only the PDF of 
one of the several fitted variables available. In principle, 
for example, we could fit the fluxes from both energy 



11 





Volume £thrcsh 


Flux N„ Nr 


a/N 




big 


30keV 


12899 345 256 


0.19% 


'Be 


std 


45keV 


9345 236 169 


0.22% 




small 


45keV 


7339 184 127 


0.24% 




big 


30keV 


42123 2256 2037 


0.16% 


PP 


std 


45keV 


27000 1335 1233 


0.19% 




small 


45keV 


21163 1025 963 


0.21% 



Table 2 

Flux en'or due to volume and threshold cuts, assuming event ob- 
served over a period 5 years. The definitions of Na, N,- and volume 
sizes are given in text. 



Volume 


^thresh 


EBe/Msc 


Epp /^pp 


Ey/Ny 


big 


30keV 


2.48% 


1.01% 


1.00% 


std 


45keV 


2.48% 


1.06% 


1.08% 


small 


45keV 


2.56% 


1.11% 


1.10% 



Table 3 

Error inherent to the likelihood method using PDFs. Es are obtained 
from the error matiix V. Ev and Ny are the error and total neutrino 
flux if the pp and 'Be are not fitted separately. 5 years of observation 
is assuined. 

and l{x) value PDFs, which, when combined, could 
lead to smaller fitting errors. However, there are signifi- 
cant correlations in fitted fluxes from energy and L (x), 
thus the combined fitting does not yield significant im- 
provement over fitting using either variable separately. 
For simplicity, we currently choose to fit using only the 
energy PDFs, which have slightly smaller errors than 
those from the L (x). 

Error induced by absolute energy scale uncertainty 
is a consequence of any energy threshold cut, because 
for events near the energy threshold, an error in en- 
ergy scale may change its status of being accepted or 
rejected. This type of error is commonly described by 
the quantity which, when multiplied by the en- 

ergy scale uncertainty 6E /E, gives the relative error in 
flux <I>. Since the absolute energy scale uncertainty is a 
property of the physical detector system, it will be ob- 
tained by insertion of a calibration source. For purposes 
of this analysis it is assumed that a 2% uncertainty on 
the scale is achievable leading to a contribution to the 
pp flux error of 0.94% at a 45keV threshold. Based on 
our simulated fluxes, values for this quantity are listed 
in Table. 4. Note that this error depends on the en- 
ergy threshold value but not the fiducial cut parameters, 
because the neutrino events are uniformly distributed 
throughout the detector volume. The calibration source 
would also be needed to establish the accuracy of the 
position reconstruction. 

These three types of errors are all considered as sys- 





'Be 


PP 


^thre.sh 


30 45 


30 45 


d*/* 
dE/E 


0.056 0.073 


0.29 0.47 



Table 4 

Coefficient of error from absolute energy scale uncertainty. 





Statistical 


Systematic 


Oscale 




'Be pp 'Be+pp 


'Be pp 'Be+pp 


'Be pp 


non-opti 


1.03% 0.61% 0.52% 


2.49% 1.08% 1.10% 


0.15% 0.94% 


optical 


1.03% 0.60% 0.52% 


2.78% 1.21% 1.19% 


~0.2% ~ 1% 



Table 5 

Relative statistical and systematic error of fitted neutrino fluxes. 
Cscalc column lists the error from absolute energy scale uncertainty 
of 2%. The terms "opticaf and "non-optical" are defined in footnote 
2. 

tematic errors, and are independent of each other, thus 
will be summed in quadrature. Table 5 lists the com- 
bined error and the statistical fluctuations of the fitted 
fluxes (without the absolute energy scale uncertainty 
part, which should add another ~ 1% to pp and com- 
bined neutrino fluxes if we assume 2% absolute energy 
scale uncertainty). 5 years worth of observation at "std" 
fiducial cut and a 45 keV energy threshold are used in 
that table. For comparison, we also show in Table 5 the 
resulting errors to be expected when all optical effects 
are included. The differences are not significant as we 
have discussed in footnote 2. 

6. Summary and Conclusions 

We have described that, on the basis of the new 
knowledge gained in recent years of neutrino properties 
and of higher energy solar neutrino fluxes, there are ex- 
cellent reasons to perform precision real-time measure- 
ments of the very low-energy neutrino fluxes from the 
Sun. The physics goals outlined in Sec. 2 include deter- 
mining the luminosity of the Sun in neutrinos, providing 
checks on some details of the SSM, testing the MSW 
effect in the LMA solution and improving constraints 
on the neutrinos mass-mixing parameters as well as pro- 
viding discovery opportunities in the new low energy 
regime. 

To achieve these goals detectors are required which 
can measure the Opp flux with a precision better than 
3% and the <I>Be flux to better than 5%. Such detec- 
tors must be capable of collecting very large event sam- 
ples and maintain good control of systematic errors. 
We have described the design of such a possible detec- 
tor, HERON, and have simulated its performance. Al- 
though the HERON detector is not presently scheduled 
or funded for construction, by experimentation in pro- 
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totypes of several liters we have measured the details 
of energy loss processes (scintillation, phonons/rotons, 
electron bubbles) for low energy electrons in the super- 
fluid [34,35,38,40] . The development of wafer calorime- 
ters capable of detecting all three channels has been 
carried out [41,42,43,44]. We have not tested an array 
of wafers as a coded aperture; however, given the well 
tested use of the method in other fields, performance as 
simulated can be reasonably expected. The superfluid 
helium target material is itself free of intrinsic inter- 
nal background and provides two channels (scintillation 
and drifted electrons) to distinguish and separate exter- 
nally entering background from point-like neutrino ES 
signals via an active coded-aperture array. 

The simulation has been directed towards establish- 
ing the systematic and other errors to be expected for 
the HERON detector in an exposure of 5 years. For that 
purpose large samples of both signal and background 
events were generated and then fully reconstructed ac- 
cording to the physical processes in helium, the de- 
tector geometry and the properties of the coded aper- 
ture design. The expected signals were based on current 
best understanding of the neutrino mass-mixing param- 
eters and the well-known electroweak scattering cross- 
sections. The backgrounds (gamma-ray conversions in 
the helium) were simulated from radioactive sources 
distributed throughout the major materials surrounding 
the helium. The level and nature of the activities as- 
sumed was in line with current best practice in solar and 
double beta decay neutrino experimentation. To sepa- 
rate a combined sample of pp, Be7 and backgrounds into 
their respective flux components, an extended loglikeli- 
hood method was used employing probability distribu- 
tion functions constructed from various samples of the 
above simulation prescription. By design, the method 
included all correlations among variables imposed by 
event reconstruction, various cuts on the data as well 
as those arising from the properties of the PDF's them- 
selves. 

The results are quite promising as can be seen in Ta- 
bles 2-5. It appears that should the detector be built and 
perform as modeled it would be capable of satisfying the 
criteria necessary for the precision pp and Be7 flux mea- 
surements. For example, to take the particular choice of 
the so-called "standard threshold and fiducial cut" on a 
5 -year exposure and combining all errors except energy 
scale (statistical, high-level cuts, likelihood method) a 
precision on pp flux of ±1.35% (or ±1.68% including 
energy scale uncertainty) results. Similarly for <I>Be, we 
find errors of ±2.96% and ±2.97% without and with, 
respectively, the energy scale uncertainty included. The 
full neutrino flux obtained without attempting to em- 



ploy the separation of <J>pp and <J>Be individually would 
present a combined error of ±1.31%. 

The validity of the background model used in these 
simulations is a key issue. The composition and rel- 
ative magnitude of the background sources assumed 
were based on current experience in the field and should 
therefore be realizable in practice; nonetheless, it is im- 
portant to ask to what extent are the simulation results 
dependent upon the assumed model and to what ex- 
tent can the model be checked in practice. The decay 
modes of the sources, branching ratios and energies of 
the decay products are well known and the method of 
their propagation through the simulation programs are 
well established. Perhaps most important is the question 
of whether we may have mis-estimated the total back- 
ground level and if so how strongly the result would 
be affected. A mistake in the magnitude would enter 
principally through the separation procedure using the 
PDF's; consequently we have tested this effect by vary- 
ing the assumed size of the background over a wide 
range. The effect is not drastic; for example, should the 
background be 5 times larger, the Opp error would dou- 
ble while a factor 50 larger background would raise the 
Opp error by six times (however, this latter rate would 
introduce prohibitive wafer deadtime). In contrast, re- 
ducing the background by a factor of 0.5 only improves 
the 4>pp error by 15%. 

In practice, there are some checks available on the 
model. Due to the good position and energy resolution 
for the point-like signal events, fiducial volumes of var- 
ious sizes can be made and the stability of the flux re- 
sults checked. As we have seen in Tables 2, 3 and 5, 
varying the high-level cuts by setting different fiducial 
volumes and thresholds does not have a strong effect 
upon the expected errors. Similarly, the nature and de- 
pendence of the observed spatial and energy distribu- 
tions as a function of these cuts can also be compared 
directly to the model. 

In conclusion we believe that a detector of the 
HERON design utilizing superfluid helium and a coded 
aperture array could provide the capability to carry out 
the multiple physics goals achievable through precision, 
real-time, simultaneous measurements of «t>pp and Obc- 
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Fig. A. 1. A typical fitting for Earth orbit eccentricity. (Inset: Nor- 
malized eccentricity results from 2 x ID'* trial simulations.) 
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Appendix A. Solar orbit eccentricity. 

In order to examine the effectiveness of our flux fit- 
ting method under a more realistic context, we chose 
to test for the annual solar neutrino flux oscillation due 
to the Earth orbit eccentricity. The orbit of the Earth 
around the Sun has an eccentricity of ~ 1 .67%. Since the 
diameter of the Sun (~ 1.4 x lO^km) is much smaller 
than the radius of the orbit (~ 1.5 x 10** km), the Sun 
can be treated as a point source, thus the neutrino flux 
observed on Earth will oscillate according to l/r^. 

We simulated the number of events observed daily 
over a span of five years. These event numbers consist 
of both pp and ^Be events including their statistical 
fluctuations. A random error according to the systematic 
uncertainty of the flux fitting method is then added to 
each day's flux, to simulate the errors introduced during 
the reconstruction and flux separation process. Then the 
daily fitted event counts are grouped into consecutive 
60 day periods, with 5 days' worth of data each year 
discarded for simplicity. Thus over 5 years, each of these 
periods contain 300 days' worth of events. A fit of 
these event numbers against the model of an elliptic 
Earth orbit can then be performed, using the eccentricity 
as the fitting parameter. 

Figure A.l shows a typical set of data for such fit- 
ting. The error bars on the "fitted" data include both sta- 
tistical and systematic errors, as discussed in Sec. 5.4, 
where the absolute energy scale uncertainty is taken as 
2%. Best fit of the particular data set in that figure gives 
an eccentricity of 2.1%; repeating this simulation for 



2 X 10"* times reveals the distribution of the best fitted 
eccentricity values to be 1.68% ± 0.95%(lo). This ex- 
ercise demonstrates that our flux fitting method is capa- 
ble of resolving the solar neutrino flux to the precision 
of a few percent given 5 years of detector running time. 
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